arXiv:1502.01800vl [cond-mat.str-el] 6 Feb 2015 


When is the Kadowaki-Woods ratio universal? 


D. C. Cavanagh, A. C. Jacko, and B. J. Powell 
School of Mathematics and Physics, The University of Queensland, QLD 4072, Australia 

We calculate the Kadowaki-Woods ratio (KWR) in Fermi liquids with arbitrary band structures. 
We find that, contrary to the single band case, the ratio is not generally independent of the effects of 
electronic correlations (universal). This is very surprising given the experimental findings of a near 
universal KWR in many multiband metals. We identify a limit where the universality of the ratio, 
which has been observed experimentally in many strongly correlated electron systems, is recovered. 
We discuss the KWR in Dirac semimetals in two and three dimensions. In the two-dimensional 
case we also generalize the KWR to account for the logarithmic factor in the self-energy. In both 
cases we find that the KWR is independent of correlations, but strongly dependent on the doping 
of the system: for massless fermions the KWR is proportional to the inverse square of the carrier 
density, whereas the KWR for systems with massive quasiparticles is proportional to the inverse of 
the carrier density. 


I. INTRODUCTION 


Fermi liquid theory describes the low temperature behavior of the vast majority of metals extremely well [IHi. 
One of the beauties of Fermi liquid theory is that it reduces the description of the interacting electron fluid to a small 
number of (Landau) parameters. Therefore, ratios in which these parameters cancel, such as the Wilson-Sommerfeld 
ratio and Wiedemann-Franz law 0,0 provided important tests of Fermi liquid theory. 

In a Fermi liquid the electronic contributions to the resistivity [pei (T) = AT'^] and heat capacity [Cei {T) = yT] 
are both governed by the effective mass, m* - roughly speaking A oc m*^ and 7 oc m*. So the Kadowaki-Woods ratio, 
Aj^"^, should be constant in a Fermi liquid [Mil- More precisely one might expect correlations to leave the Kadowaki- 
Woods ratio (KWR) unrenormalized because a Kramers-Kronig transformation relates the real and imaginary parts of 
the self-energy |lll - ll3l| . which determine the electronic contributions to the heat capacity and resistivity respectively. 
This means that the KWR is somewhat similar to a fluctuation-dissipation theorem. 

First Rice (a and later Kadowaki and Woods [l^ found that Ajj^ is approximately constant within classes of 
materials (transition metals and heavy fermion compounds, respectively). However, the ratio differs by two orders 
of magnitude between these two classes. It was subsequently discovered that the Kadowaki-Woods ratio (KWR) in 
transition metals and organic charge transfer salts can be even larger than in the heavy fermions (see Refs. [Hill 
and references therein.) 

It was l ong believed [ll|, [13 that the size of the KWR gave an indication of the strength of the electron-electron 
scattering ^ |. However, this has been shown to be incorrect [a [ 3 . Rather, the large variations in the KWR between 
different classes of materials can be explained almost ent irely by taking into account non-interacting properties of the 
materials (e.g. electron density and dimensionality) [l^ [i0- Furthermore, it has been shown [13 that the modified 
KWR A/y^/, where / is a material specific function of the non-interacting band structure (defined below), takes 
the same predicted value, SI/dTr/ifcge^, in a large range of transition metals, charge transfer salts, heavy fermion 
compounds, and elemental metals, a result which has since been verified in many other materials | 18 l - l 22 j| . 

Deviations from this universal predicted value of the KWR could provide an indication of non-Fermi liquid behavior. 
Most previous calculations of the (modified) KWR BEO have focused on simple, single-band models with toy 
dispersion relations, e.g. spherical Fermi surfaces. Understanding the effects of more complicated and realistic band 
structures is important if the modified KWR is to be used to identify deviations from Fermi liquid theory. 

Previous studies of the KWR in systems with orbital degeneracy [ 2 ^, or multiple bands [l3 have found that 
either of these can cause significant variation in the KWR. In this paper, we derive a modifed KWR for Fermi liquids 
with arbitrary dispersion relations, including band structures with multiple bands. We find that the universality of 
the KWR evident in the single band expression is not a general feature of the multiple band case. In particular in the 
most general case the strength of electronic correlations does affect the value of the KWR. This is extremely puzzling 
as the KWR is found to be close to its universal value in many multiband systems. However, if the renormalization 
is the same on all bands (in a sense made precise below) correlations do cancel from the KWR. 

In the case of uniform renormalization across all bands, we demonstrate, that the KWR Rkw oc l/niV^ for massive 
quasiparticles in a system of N}, bands. This is particularly relevant to semimetals where the low career density opens 
the possibility of large variations in the career density, n. This is further enhanced in Dirac semimetals - we find that 
the massless fermion dispersion relations lead to Rkw oc l/ri^N^ in three dimensions. In two dimensions a similar 
result holds once the KWR is generalized to account for the logarithmic factor in the imaginary part of the self-energy. 
Furthermore, we show that the logarithmic factor in the imaginary part of the self-energy leads to an increase in the 
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KWR, which may provide an experimental signature of this factor in the self energy. 

The remainder of this paper is laid out as follows. In the following section, we calculate the resistivity of an arbitrary 
multi-band Fermi liquid by calculating the conductivity from the Kubo formula of linear response. In Section fllll we 
calculate the effect of multiple bands on the heat capacity. In Section El we combine these results to determine the 
form of the Kadowaki-Woods ratio in arbitrary band structuress. In Section |Vl we then apply this form for the KWR 
to simple models of Dirac semimetals in two and three dimensions. 


II. CONDUCTIVITY FROM THE KUBO FORMULA 


In general, the contribution to the intraband self-energy from interband terms scales quadratically with the intra¬ 
band self-energy, such that, as long as the self-energy is small, the intraband contribution will dominate [1^. When 
intraband scattering is the dominant contribution to the scattering rate, the diagonal component of the conductivity 
tensor, axx, for a material with Nb bands crossing the Fermi surface in the low temperature limit is 
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where Saa (w) = S'aa (w) -t-iS-'a (w) is the self-energy, = \l- (w) /duj]ui=o\~^ is the quasiparticle weight, vuxa 

the x-component of the group velocity of an electron in spin-band d = (a, a), where a denotes the band and a the 
spin (giving 2Nb spin-bands), and we have neglected vertex corrections (cf. Ref. El). The sharply peaked derivative 
of the Fermi-Dirac distribution at low temperatures implies that 
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where (...) indicates an average over the Fermi surface. The conductivity is clearly then the sum, in series, of the 
conductivities of the individual bands: 
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where the bare density of states of band d at the Fermi level is Dg-a = (27r)“^ d^k [d (wka — m)] <^xxa is the 
conductivity of due to spin-band d. It immediately follows that, if interband scattering is neglected, the resistivities 
of the bands must add in parallel, consistent with Matthiessen’s rule [I|. 


1. The Quadratic Contribution to the Resistivity 


To determine the explicit form of the conductivity in a Fermi liquid, we generalise a local (momentum-independent) 
phenomenological self-energy model proposed by Miyake, Matsuura and Varma, m to include multiple Fermi surfaces. 
We assume that in each band the self-energy takes the phenomenological form 


S"a(^) = - 
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(4b) 


where the energy oo is measured from the Fermi level, is the impurity scattering rate, sa = ''^a/37rZ?o;a is the 

electron-electron scattering rate in the unitary limit (with na the density of charge carriers in band d), F{x) 

is a monotonically decreasing function with F"(l) = 1 and F{oo) = 0 and w- is an energy scale characterising the 
strength of the many-body correlations (we will show below that Z^ ~ a;t/4sa). This is the natural form for a 
local Fermi liquid provided interband interactions are weak compared to the intraband interactions. If the two are 
comparable (as in some multiorbital models) the self energy is proportional to Nb — 1 23|. 
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In the limit of vanishing impurity scattering, Tq,^ 0, the conductivity is then given by 
^ ^ Dl.a{vlxa) f 1 dnfjuj) 

^2 na J ^w2 + (7rfcsT)^ duj 


After computing the energy integral, the A^^-band A coefficient is 

\2 \ 
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thus, we see that the coefficients of the individual bands, Aa, add in parallel. 


(5) 


( 6 ) 


III. THE HEAT CAPACITY VIA KRAMERS-KRONIG TRANSFORM 


It follows from the extensivity of the heat capacity that, in a multi-band system, the total heat capacity is given 
by the sum in series of the heat capacity due to each individual band [I|, 
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To determine how the heat capacity is influenced by interactions, and calculate the relevant 7 coefficients, we 
need to first find the real part of the self-energy in each band. The Fermi liquid self-energy, being a causal res pon se 
function in the time domain, satisfies the conditions for the Kramers-Kronig relations in the frequency domain [l2| - 
in particular lim S" (w) —>■ 0. Knowledge of the form of the imaginary part of the self-energy is therefore sufficient 

Ld—¥CO 

to determine the real part, which appears in the definition of the quasiparticle weight and therefore in the expression 
for the heat capacity. The real part of the self-energy within each band is then 


00 



where we have again taken the limit of vanishing impurity scattering, and have restricted the pole of the integral, 
w' = uj, to occur below the cut-off energy scale, \u:'\ < |w|| (i.e., only considered low-energy excitations). The first 
term in Eq. (|S]) contributes a logarithmic term to the result, which we approximate by the lowest order terms in a 
Taylor series exjransion, while the second and third terms contribute linearly to the self-energy. Neglecting terms of 
order O (w/w-) and higher, we find that the real part of the self-energy for a low-energy quasiparticle in band a is 


where 


S'aaH = - 


driaw 






(9) 


( 10 ) 


it follows straightforwardly from the definition of F{y) that 1/2 < ^ < 1, henceforth we take ^ = 1 for simplicity. 
Inserting this expression for the real parts of the self-energies into the heat capacity expression Eq. o, and taking 
the strong scattering ( 4 spa/ 7 rw| 1 ) limit, which corresponds physically to m* 2 > toq, we obtain 
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Note that Eqs. © and (El) imply that Wj ss 4sa-^a = 4Zana/37ril)o;a, which gives a straightforward interpretation of 
this energy scale. 
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IV. THE KADOWAKI-WOODS RATIO 


Given the above calculations of A and 7 for the multiple band system, Eqs. 
is given by 


and (ED, we find that the KWR 
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Alternatively, one may write 
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where we have defined the material specific function for an V;,-band Fermi liquid with the resistivity measured in the 
X direction, 


/x;Ar, ({<} , {na}) = 


r2Nb 1 

2 



_d—l 



2 Nt, £>•2 


i (7)‘ 


6=1 


(14) 


In the one-band limit, this expression simplifies to that calculated in Ref. [13 (with = nj, = n/2 and = 
-^0,4, = -Do/ 2 ) 


/x;i (n) — n{vQ^)DQ. 


(15) 


The inclusion of multiple Fermi surface sheets significantly complicates the form of ({w^}, {na}). Most 
importantly which describes the electronic correlations, does not cancel out of the multiband expression as it does 
for the single band KWR [l^. Therefore, our calculation predicts that the Kadowaki-Woods ratio is not, in general, 
independent of electronic correlations. This is rather surprising as observed values of the KWR (including the values 
for many multiband systems) are in almost universal a gree ment with the prediction from the single band calculation, 
that electronic correlations do not influence the KWR |1^ . 

It is therefore important to ask how renormalization effects might cancel in the multiband case and hence universality 
might be recovered. 

The simplest limiting case for which the effects of many-body correlations cancel out of the KWR is when is 
independent of the band index, a. This is a straightforward extension of the earlier assumption of the locality of the 
self-energy, by assuming that it is independent of band index as well as momentum. This assumption yields 
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Though this calculation has been performed with exactly uniform correlation strengths for simplicity, the result will 
hold approximately while the correlation strengths are close to uniform. 

Other limits do produce a universal KWR, for example, for free fermions in three dimensions Dp.g oc then if 
the carrier density in one band is much larger than all others the correlations cancel from the KWR. But it seems 
unlikely that this is relevant the behaviour of a broad range range of materials. Note, in particular, that if we have 
a single heavy band it will dominate the heat capacity, but be shorted out of the resistivity. Therefore, the limit of 
a single heavy band is far from universal. Therefore, the above calculation seems to suggest the correlation strength 
(as measured by w- = iZ^nd/^T^Do-d) does not vary strongly between different bands in strongly correlated systems. 

If, further, all of the bands are identical, i.e., if the carrier density n, the Fermi velocity and density of states 

Dq are equal for all bands, we have 


fx-Nb (n) = Nl{vl^)nDl = N^fdx-,i{n), 


(17) 


and 


n _ 81 Rkw-,1 


( 18 ) 
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Thus we see that the expression for the KWR in the single-band case is modified by a simple factor of 

At first glance this expression appears rather similar to the finding of Kontani et al. [231124| that in the multiorbital 
periodic Anderson model with No impurity states the KWR is reduced by a factor of l/No{No — 1). However, 
closer examination reveals that the results are actually very different. In particular Kontani’s factor arises because 
of interorbital terms in the self energy whereas our 1/A^^ factor arises purely from the electronic structure. We do 
not obtain Kontani’s factor because of our assumption that interband interactions are irrelevant at low energies. In 
contrast the model Hamiltonian studied by Kontani explicitly sets the intra- and inter-orbital interactions to the 
same strength. Which approach is appropriate will depend on the material. This therefore adds another layer of 
non-universality to the KWR. 


V. DIRAC SEMIMETALS 

The general expression, Eq. (USD, can be applied to systems of arbitrary band structure to calculate the generalised 
Kadowaki-Woods ratio, taking into account the effects of multiple bands. Even for materials where the quasiparticle 
weight is similar for all bands such efforts will, in general, involve first principles band structure calculations. In this 
section, we calculate the KWR for a simple, linear dispersion (eta = |k|) appropriate for Dirac semimetals. These 

models present analytically tractable and instructive examples of complicated band structures for which the presence 
of multiple bands is important (here, the number of bands may be taken as equivalent to the number of Dirac cones, 
as each cone will form a sheet in the Fermi surface). Furthermore, these materials are an example of materials where 
the different sheets of the Fermi surface where we expect the correlations to affect all bands equally and the low carrier 
density suggests that correlations will play an important role. We also determine an expression for the KWR in two 
dimensions, accounting for the logarithmic factor arising in the self-energy [25l - [^ . 


A. 3D Dirac Semimetals 

We first apply the KWR expression to a simple model of a two-band three-dimensional Dirac semimetal, e.g. 
Cd 3 As 2 , which possesses properties similar to doped graphene The spatial symmetry between the two bands 

and the spin symmetry simplifies the calculation greatly; applying the expression for Ni, identical bands, Eq. (ED, 
we find 

7^ f3x,2 47re2te| ’ 

where f 3 x ,2 = ‘In?/-Kk? is the material specific function for a 3D Dirac semimetal and is larger by a factor of 4 than 
the corresponding / in a single-band calculation with the same dispersion. We note that the Kadowaki-Woods ratio 
for these materials depends straightforwardly on the electronic density, which is tunable via chemical doping , 

providing a potential experimental test of this expression. 


B. Doped Graphene and the Self-Energy in Two Dimensions 


In calculating Eq. ED, we utilised the form of self-energy given in Eq. (SD, which implicitly assumes a three- 
dimensional material. In two dimensions, the self-energy differs from the three dimensional case, with the inclusion 
of an additional logarithmic factor [i,|3jli3. In order to calculate the KWR for doped graphene, we introduce the 
following model for the imaginary part of the self-energy at low energies and temperatures for 2D systems 


= 



( 20 ) 


where Ba is a constant of order unity (for example the calculation in Ref. |4l| gives Ha ~ I/tJ" for graphene; but other 
calculations give slightly different results). We stress that, in the zero and high temperature (fcsT » w) limits, Eq. 
(ED reproduces the known results for those limits [13, li^|4l|. We further assume that, above the relevant energy 
scale, w*, the self-energy decreases monotonically, as in the 3D case. From this expression for imaginary part of the 
self-energy, we calculate the real part using the Kramers-Kronig transformation, and find that 
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where the logarithmic factor arises due to the requirement that the self-energy be continuous at w = w*. 

In defining the 2D KWR, we must account for the fact that the logarithmic contribution to the imaginary 
part of the self-energy (and scattering rate) results in a corresponding logarithmic factor in the resistivity, p = 
—AT^ log {ttUbT/B uj*). The coefficient of this logarithmically adjusted quadratic term will be used in our expression. 
Calculating the conductivity from the Kubo formula gives 


A. = 


8kg 


T^e^hna{vl-^)ujlDlB 


( 22 ) 


where, in evaluating the energy integral (see Eq. ([5])), we have approximated the logarithmic contribution by its value 
at the Fermi surface. From the derivative of the real part of the self-energy, we find the value for the quasiparticle 
weight and therefore the linear coefficient of the heat capacity: 


la 


27rfc|na 
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Defining the two-dimensional Kadowaki-Woods ratio as we find 


Rkw = —^ = 
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is the two-dimensional material specific function, which for graphene takes the form 


f 2 x -2 (n) = 4/2 x;i (n) 


2n? 

irh? 


log (B) + i 


(25) 


Here the Kadowaki-Woods ratio again depends on the electron density in a straightforward manner, but differs from 
the previously derived expression (Eq. ([19])) only by a factor of [log (B) + 1/2]^. Assuming B ~ I/tt, this term more 
than doubles the expected Kadowaki-Woods ratio. Measuring the KWR of graphene presents a number of technical 
challenges, chiefly the difficulty of performing calorimetric measurements on a material of single atom thickness. It 
has been a long standing problem to observe the logarithmic factor in the resistivity of a two-dimensional Fermi liquid 
because it would require data over many orders of magnitude in temperature. However, observation of this correction 
offers a potential test of the dimensionality of the Fermi liquid self energy without the need to take data over multiple 
decades. 


VI. CONCLUSIONS 

We have shown that, in general, the Kadowaki-Woods ratio of a multiband local Fermi liquid is changed by electronic 
correlations. This is in marked contrast to the single band case, where the KWR is independent of the strength of the 
electronic correlations. It is therefore puzzling that the experimental data suggest the within classes of materials the 
KWR is remarkably consistent, and that the modified KWR is remarkably consistent across many chemically diverse 
strongly correlated metals. The simplest explanation is that the correlations are indeed very similar across all bands 
in these materials. We have also shown that a non-parabolic dispersion does not significantly alter the form of the 
KWR, provided the fermions remain massive. 

In the case of uniform renormalization across bands, we have further demonstrated that Rkw ot 1/nN^ for massive 
quasiparticles in a system of Nb bands. This is particularly interesting in semimetals where the low career density 
opens the possibility of large variations in the career density, n. This is further enhanced in Dirac semimetals, where 
we have shown that the massless fermion dispersion relations lead to Rkw oc 1/n'^N^ in three dimensions. In two 
dimensions a similar result holds once the KWR is generalized to account for the logarithmic factor in the imaginary 
part of the self-energy. Furthermore, we have shown that the logarithmic factor in the imaginary part of the self-energy 
leads to an increase in the KWR, which may provide an experimental signature of this factor in the self energy. 
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When is the Kadowaki-Woods ratio universal?: Supplementary Material 


I. SCATTERING AND THE SELF-ENERGY IN ARBITRARY BANDSTRUCTURES SYSTEMS 


We wish to determine the form of the self-energy in arbitrary bandstructures and the effect of multiple Fermi 
surfaces on the self-energy and Green’s function. The elements of the dressed Nb x N}, Green’s function matrix are 
then defined by the matrix form of Dyson’s equation [l| and we have (here and below, we use the convention b ^ a 
to distinguish intraband and interband contributions) 





^a,b (^) • 



The inverse can be calculated using the adjugate matrix ii (A^ (adj [A]) ), each element of which contains 

JVb — 1 terms, with a leading term of the form Y\ \ (w), with all other terms of order O ( ^ (w) J and higher. So, 

with minimal interband scattering, the intraband Green’s function elements will be dominated by the intraband self¬ 
energy contributions, and the individual bands will reduce to separable channels. We consider in detail the two-band 
case below, to give an indication of how this separability arises.]^ 


A. The Two-Band Case 


The inverse of the two-band Green’s function matrix has the simple representation 
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G{uj) = 


[G(a.)] 


-1 


G 2 ^ ^1.2 (w) 

E2.i(cc) Gr'(cc) 


( 1 ) 


where G^j ^ (w) = Gi°a (w) — '^a,a (w) is the dressed Green’s function in band a neglecting the effects of interband 

scattering (i.e. 'Ea.b (w) = 0), and the determinant is 


[G (W)]-^ = n [Ga,a (W)]-^ - Si .2 (CC) E 2 .I (uj) . 


( 2 ) 


The intraband elements of the dressed Green’s function matrix are then given by 


Ga,a (w) — 


G-' (w) - 


(^) ^b,a (^) 

gTm 


-1 


(3) 


and the interband elements by 


Ga,b (w) — 


n G-i (w) 

n—1 

^a,b (^) 


^5,a (^) 


-1 


(4) 


The interband self-energy elements appear here as corrections to the single-band form of the intraband Green’s 
function elements, G^ „ (w). In the limit of small interband scattering, these corrections become negligible, being 

of order O ^|Ei ^2 As the interband elements of the self-energy approach zero, so do the Ga,b (w) elements, 

as interband propagation occurs only in the presence of interband scattering. Therefore, in the limit of vanishing 
interband scattering, the dressed Green’s function matrix is diagonal (i.e. intraband only), 

lim Ga,a (w) = Ga (w) , 

lim Ga,b (w) = 0. 
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The matrix elements Eqns. © and 0 can be expanded by expressing each self-energy as a sum of its real and 
imaginary parts, and making use of the knowledge that the self-energy matrix must be Hermitian (so (w) = 
Sf,,a (<^)^)- The diagonal component is then 


Ga,a (w) — 


(a;)l - E' , {to) - *E'' (w) 


ISa.h (w)r 


-1 




(5) 


For notational ease, we define ?7k,a (w) (the renormalised real part of the inverse Green’s function) and Tq (w) (the 
inverse quasiparticle lifetime) as 


(w) = w - ek,a - K,a > 

Ta (w) = E" „ (w) - (5, 


( 6 ) 

(7) 


where ek:a refers to the dispersion for the a band, measured from the Fermi energy. We find a more explicit 
expression for the intraband Green’s function element 


Ga,a (w) 


'rgyy+Tl (a;)] 

’Ik.a (w) 1 

Jlk,b M + Tb 

- |Sb,a (w)|^77k,b (w) 


?7k,a (w) 

^^k,b + Tfe (w)^ 

* [^k,b y) 

— Eb,a 

+Tb y) 

y)i^yk,b 

[Ta (W) 

(uj) 

2 

+ 

b(w) 

Ta H ( 
+ Tb (w) 

^k,b(‘^)+rb (w)^ 
) +|Eb,a (w)I^Fb 

+ 1 

(w) 

Cb,a (w)l'Fb (CU) 

2 

^k,a (w) 

jii.b y)+ Tb 

- |Eb,a (w)|^77k,b (w) 

2 

Ta H ( 

^k,b(‘*^) + rb 

+ |Sb,a (w)| Fb (w) 

2 


( 8 ) 


and, for the interband elements of the dressed Green’s function matrix, 
Ga,b{0jy^ = 


2 

n 

Gn}n (w) 

- K,n (w) - y) 

-|Sa,b(w)|" 

n—1 





Explicitly expanding the product over the bands, we obtain from this the expression for the interband elements of 
the Green’s function 

^a.h y) (jlKa y) »7k.b - Ta (w) Ff, (w) - \T,a,b (w)|^) - E" j, (w) (F^ (w) 77k,6 (w) -f 77k, a (w) Ff, (w)) 


Ga,b (<jj) — 


|Sa .6 (W)r 


i (^(Fa (w)77k,b (w) -I- 77k,a (w)Fb (w)) (w) -I- E" J, (w) (rik^a (w)77k,b (w) - Fa (a;)Fb (w) - |Ea,b (w)|^^^ 


|Sa,b yy 


(9) 


which, as expected, is proportional to the interband self-energy contribution, so that, in the limit of vanishing interband 
scattering, the interband Green’s function elements vanish also (i.e. lim Gab = G^^l = 0). 

t-s-n ’ 


B. The Spectral Density Function 

For the case of a two-band material, we concentrate here on the intraband (diagonal) terms Ga,a (w) (Eq. ([5])), 
taking the low energy limit [l[ , as the interband terms do not arise in the Kubo formula for conductivity (see below). 
In the more general case, for small interband scattering similar arguments to those above hold and the individual 
bands are separable, with the intraband spectral density functions reducing to their individual single-band forms. For 
the intraband elements in the two-band system, the spectral densities are 
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-^k.a (‘^) — 


-2 


(’7k,b + Tfc M) Ta (w) (?7k,b (‘^) + rg (w)) + |Sa,b (w)l^ Tb (w) 


’7k,a (w) (f7k,b (‘^) + Tfc (w)) - |Sa,b (w)l^ ?7k,b (w) + Ta (w) (t?^ (w) + (w)^ + |Sa,b (w)l^ Tb (w) 


( 10 ) 

Again, we take the low energy limit, in the two band case given by Fa^b = S"0, restricting quasiparticles in both 
bands to be close to the Fermi surface. We discard terms of order O (s"^b (‘^)) higher to obtain 


Ak,a (^) — 


-2’7k,b Ta ’7k,b + Sa,b i^)^b (w) 


( 11 ) 


’7k,a (w) ,, (w) - (w) ?7k,b (w) + F^ (w) (w) + (w) Fb (w) 

We can simplify this expression further by assuming that the interband self-energy E^^b (w) is much smaller than the 
single band self-energies, which gives the low-energy limits of the spectral density functions, 


lim Ak a (w) = lim — 
r->o ’ r-i-o '■ 


-2’7kb(‘*’)ra (‘^) 


’7k,a (w) rjl i^ (w) 
-2Ta{uj) 

= iim —7^—^^—- 

+ Ta (w) 

= -2ttS (?7k,a (w)) , 
-2TTd (ryk,a (w)) 


(‘^)’7k,b (‘^) 


( 12 ) 

(13) 


lim A,a (w) = , , 

r^o Fa(a;) 

It can be seen from Eqn. (1131) that, in the limit of small interband scattering, the spectral density function for 
band a depends only on the properties of band a and the two bands represent independent channels. In any case, in 
the conductivity derivation below, the interband terms do not play a role, as the current-current correlation function 
restricts the currents to the intraband case. 


II. INTERBAND SCATTERING 


The spectral density function for the interband (off-diagonal) elements can be expressed in the low energy limit as 

Sa,b(‘^)(ra(w)ryk,b(w)-f77k,a(w)rb(w)) K,b 


Ak,a,b (^) — 


?7k,a (w) ryk,b (w) - E^^^ ’ 


(14) 


(’7k,a (w) ?7k,b (w) - (^)) 

which clearly vanishes in both of the limits E^^b (w) —0 (no interband scattering) and Ea,b (w) —>■ oo (strong interband 
scattering). 


A. Effects of significant interband scattering on the intraband self-energy components 


If we avoid making the approximation that the interband elements of the self-energy are vanishingly small, we 
have the expression for the intraband spectral density functions Eqns. (HD, which can be rearranged to yield the 
expression, for b ^ a 


-2’7k,b (w) 


Ak,a (w) — 


F, 




vtb H 


(w) - "Ik.b (w) 


■ <,b (^) 


i,a (w) -f ( ) Fb,b (w) 


-2 


Ha(^) + (|;^)'rb.b(a.) 


2 r n 2 ’ 

VKa (w) - ( ) ’7k,b (w) -f Fa,a (w) + ( ^(a;) ) (‘*’) 


(15) 
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If we again take the low energy limit, we obtain the expressions for the spectral density function 


lim Ak.a (w) = -27r5 ( ?7k,a (w) - ( ] ) ??k,6 (w) 

r^o \ V^k.&(w) ' 


and spectral density function squared 


2ttS ("r^k.a (w) - Vk,b (w) 


lim Ak „ (w) = 

ha,a (w) + ) ^>>,b (w) 


The two band conductivity (see below for derivation of the conductivity formula) is then 



which can be seen to approach the form given in the main text for |S']^ 2 | ^ l^k.ol, and, in the opposite limit 
1 ^ 1 , 2 ! \Vk,a\ the corresponding delta function takes the form <5 2 /'' 7 k,a)^ ?7k,a^ and it is obvious that argument 

of the delta function will never vanish in this limit, and there will be no contribution to the conductivity from the 
a band. The intermediate range is more complicated, with the limit 151(^2 1 ~ l^k,a| giving an expression for the 
conductivity 


(jI'5 = he^ 


d^k 


(27r) 


3 l^k,a:| 


duj 

(27r) 


-25 (?7k4 {(Jj) - r]k ,2 (cu)) 

[hi,! (w) + r2,2 (w)] 


dnfjuj) 

dw 


so the conductivity will vanish unless the two bands’ renormalised energies overlap. 


(16) 


III. DERIVATION OF THE CONDUCTIVITY FORMULA 

In a multi-band system, the electronic current density operator component in the a direction ja, and the current- 
current correlation function iTa-y (iojn) are 


ja — e 'y ( "CkjCejaCk ,jCk,(T,a) (1"^) 

k.(7,a 

/3 

TTaj (ilMn) = ^ j {Trja [t) (0)), (18) 

0 

where, Vo is the U-dimensional volume of the unit cell, —e is the charge on an electron, cj^ ^ and Ck,a respectively 
create and annihilate an electron with momentum hk and with spin and band given by the spin-band index d = (a, a), 
yk,a,a = ^~^dek,d./dka is the a-component of the (group) velocity of an electron with momentum Hk in band a, r is 
the imaginary time (/3 = {kBT)~^), is the imaginary time ordering operator, and the nth fermion Matsubara 
frequency, with (...) denoting an ensemble average over the system. The lowest order (non-interacting) term in an 
S'-matrix expansion for the correlation function is given by 
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4°7 = 




dre ” 'yk,a,aWk', 7 ,a'( 7 TC|{j ''k,cr,a (t) (0) Ck’,cr',a' (0)). (19) 

k,k’;(T,cr' ;a,a' 

Using Wick’s theorem we then rearrange the terms in the correlation fnnction Eqn. dm), to express it in terms of 
non-interacting Green’s functions [ij 

a (r) Ck,a,a (t) ( 0 ) Ck’,^'.a' ( 0 )) = (ck,a,a (t) 4 ,( 0 )) (ck’,a',a' ( 0 ) c]^ (t)) = “T) , 

where the ensemble averages enforce the conservation of momenta, as well as spin-band index. The non-interacting 
correlation function is then 


^7 = «k.a,ank. 7 .a J cfl (k; r) cfl (k; -r), 


( 20 ) 


to which we introduce interactions by replacing the bare Green’s functions with their corresponding dressed form [l| , 
taking interactions into account solely via the self-energy (i.e. neglecting vertex corrections). We then perform the 
Fourier transform to give the Matsubara frequency correlation function (from which we can find the finite temperature 
form via analytic continuation) 

1 ^ 


g2 I 

^ct 7 ^ ^ ^k,a,a^k, 7 ,a ~q ^ ^ Gq^q (k, ip -\- icu') Ga^a (k, ip^ . 

k.a ip 


( 21 ) 

The sum over Matsubara frequencies is performed by moving to the spectral representation, where the interacting 


equenc 

Green’s functions are given by [ 11 , Ga,d. (k;ia;) = f (k, w) / {ihuj — co), with Aa the spectral density function 
for an electron in spin-band d. In the limit of weak interband scattering, the interband contributions to A^ vanish 
quadratically (see Supplementary Material for a detailed derivation), and the spectral density function is approxi¬ 
mately that of a quasiparticle in the d band, with only the intraband contributions contributing significantly. We 
then find 

S'a.a(*Wn) = g X! ^ 

ip 


where [l| 


Thus, 


duji 

2tt 


Aa{k,uii) j ^Aa{k,ui2)S^^l{iLJn), 


^ (k;.p + ™)G<"> (k;.p) = 

o — inuin + UJ 2 — UJi 


IP 


/ duj 

-^Aa (k,a; 2 ) Aa (k, a ;2 -I- w) {nf{uj2) - nf{u)2 + h^)) ■ 


( 22 ) 

(23) 

(24) 


The conductivity is then found by taking the zero frequency limit of the retarded current-current correlation function 
divided by the frequency w, 

Im {'Ka~i (w + *(5)} 


J Q,^ 


= lim 
5->0 
cj —>^0 


UJ 

CXD OO 

2tV^ f d^k f doj 2 n dnf{u})\ 

“ ® i ^^^’k.a.a«k,7.a J —Aa{k,uj) ^j . 


(25) 


The low energy limit of the spectral density function with weak interband scattering is given by lim A^ ,di^) 

Im{S}—>-0 

—2'!rZS [oj — Zei^-a) /U~ g (w), with the quasiparticle weight in band a. 


Zji = 


1 - 


d^'d.d (^) 


duj 


(26) 


|k|=feF 


The A^fe-band conductivity is then 


OO OO Q 7^ 

r d^k f duj-2TTZa6 {uj - ZaS-k-a.) 2 f dnf{uj)\ 

J ^ J 

-00 —00 
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